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Abstract 

Practically all serious calculations of exclusive particle production In ultra-relatMstlc 
nuclear or hadronic Interactions are performed In the framework of Grlbov-Regge the- 
ory or the elkonaUzed parton model scheme. It Is the purpose of this paper to point out 
serious Inconsistencies In the above-mentioned approaches. We will demonstrate that 
requiring theoretical self-consistency reduces the freedom in modeling high energy 
nuclear scattering enormously. 

We will Introduce a fully self-consistent formulation of the multiple-scattering scheme 
In the framework of a Grlbov-Regge type effective theory. In addition, we develop new 
computational techniques which allow for the first time a satisfactory solution of the 
problem in the sense that calculations of observable quantities can be done strictly 
within a self-consistent formalism. 



1 Open Problems 

With the start of the RHIC program to investigate nucleus-nucleus collisions at very high 
energies, there is an increasing need of computational tools in order to provide a clear 
interpretation of the data. The situation is not satisfactory in the sense that there exists 
a nice theory (QCD) but we are not able to treat nuclear collisions strictly within this 
framework, and on the other hand there are simple models, which can be applied easily 
but which have no solid theoretical basis. A good compromise is provided by effective 
theories, which are not derived from first principles, but which are nevertheless self- 
consistent and calculable. A candidate seems to be the Grlbov-Regge approach, and 
- being formally quite similar - the eikonalized parton model. Here, however, some 
inconsistencies occur, which we are going to discuss in the following, before we provide 
a solution to the problem. 

Grlbov-Regge theory |[l|, |2| is by construction a multiple scattering theory. The ele- 
mentary interactions are realized by complex objects called "Pomerons", who's precise 
nature is not known, and which are therefore simply parameterized, with a couple of 



1 



parameters to be determined by experiment ||^] . Even in hadron-hadron scattering, sev- 
eral of these Pomerons are exchanged in parallel (the cross section for exchanging a 
given number of Pomerons is called"topological cross section"). Simple formulas can be 
derived for the (topological) cross sections, expressed in terms of the Pomeron parame- 
ters. 

In order to calculate exclusive particle production, one needs to know how to share 
the energy between the individual elementary interactions in case of multiple scatter- 
ing. We do not want to discuss the different recipes used to do the energy sharing (in 
particular in Monte Carlo applications). The point is, whatever procedure is used, this 
is not taken into account in the calculation of cross sections discussed above So, 
actually, one is using two different models for cross section calculations and for treating 
particle production. Taking energy conservation into account in exactly the same way 
will modify the (topological) cross section results considerably. 

Another very unpleasant and unsatisfactory feature of most "recipes" for particle 
production is the fact, that the second Pomeron and the subsequent ones are treated 
differently than the first one, although in the above-mentioned formula for the cross 
section all Pomerons are considered to be identical. 

Being another popular approach, the parton model amounts to presenting the 
partons of projectile and target by momentum distribution functions, fi and fj, and 
calculating inclusive cross sections for the production of parton Jets as a convolution of 
these distribution functions with the elementary parton-parton cross section d(Tij/dp\, 
where t, j represent parton flavors. This simple factorization formula is the result of can- 
celations of complicated diagrams and hides therefore the complicated multiple scatter- 
ing structure of the reaction, which is finally recovered via some unitarization procedure. 
The latter one makes the approach formally equivalent to the Gribov-Regge one and one 
therefore encounters the same conceptual problems (see above). 

2 A Solution: Parton-based Gribov-Regge Theory 

As a solution of the above-mentioned problems, we present a new approach which we 
call "Parton-based Gribov-Regge Theory": we have a consistent treatment for calculat- 
ing (topological) cross sections and particle production considering energy conservation 
in both cases; in addition, we introduce hard processes in a natural way. 

The basic guideline of our approach is theoretical consistency. We cannot derive ev- 
erything from first principles, but we use rigorously the language of field theory to make 
sure not to violate basic laws of physics, which is easily done in more phenomenological 
treatments (see discussion above). 

Let us first introduce some conventions. We denote elastic two body scattering am- 
plitudes as T2^2 and inelastic amplitudes corresponding to the production of some final 
state X as T-^^x (see fig. As a direct consequence of unitarity on may write the 

X 

Figure 1: An elastic scattering amplitude T2^2 (left) and an Inelastic amplitude T2_,x (right), 
optical theorem 2Imr2^2 = Yl,x(^'^^x){T2^x)* ■ The right hand side of this equation may 
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be literally presented as a "cut diagram", where the diagram on one side of the cut is 
(T2^x) and on the other side (T2^x)*. as shown in fig. ^. So the term "cut diagram" 




Fignre 2: The expression y^ y(T2-^jy).(T2^x)* which may be represented as a "cut diagram". 

means nothing but the square of an inelastic amplitude, summed over all final states, 
which is equal to twice the imaginary part of the elastic amplitude. 

Before coming to nuclear collisions, we need to discuss somewhat the structure of 
the nucleon, which may be studied in deep inelastic scattering - so essentially the 
scattering of a virtual photon off a nucleon, see fig. ^. Let us assume a high virtuality 




Figure 3: Deep Inelastic scattering: the general diagram (left) the "sea" contribution with a soft Pomeron at 
the lower end (middle) and the "valence" contribution (right) . 

photon. It is known that this photon couples to a high virtuality quark, which is emitted 
from a parton with a smaller virtuality, the latter on again being emitted from a parton 
with lower virtuality, and so on. We have a sequence of partons with lower and lower 
virtualities, the closer one gets to the proton. At some stage, some "soft scale" scale 
Qo must be reached, beyond which perturbative calculations are no longer valid. So 
we have some "unknown object" - indicated by "?" in fig. |^ - between the first parton 
and the nucleon. In order to proceed, one may estimate the squared mass of this 
"unknown object", and one obtains doing simple kinematics the value Qf^/x, where x is 
the momentum fraction of the first parton relative to the nucleon. Therefore, sea quarks 
or gluons, having typically small x, lead to large mass objects - which we identify with 
soft Pomerons, whereas valence quarks lead to small mass objects, which we simply 
ignore. So we have two contributions, as shown in fig. ^: a sea contribution, where the 
sea quark or gluon is emitted from a soft Pomeron, and a valence contribution, where 
the valence quark is one of the three quarks of the nucleon. The precise microscopic 
structure of the soft Pomeron not being known, it is parameterized in the usual way a a 
Regge pole. 

Elementary nucleon-nucleon scattering can now be considered as a straightforward 
generalization of photon-nucleon scattering: one has a hard parton-parton scattering 
in the middle, and parton evolutions in both directions towards the nucleons. We have 
a hard contribution Thaid when the the first partons on both sides are valence quarks, 
a semi-hard contribution Tgcmi when at least on one side there is a sea quark (being 
emitted from a soft Pomeron), and finally we have a soft contribution, when there is no 
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Figure 4: The soft elastic scattering amplitude T^oft (left), the hard elastic scattering amplitude Thard (middle) 
and one of the three contributions to the semi-hard elastic scattering amplitude Tgemi (right) . 

hard scattering at all (see fig. We have a smooth transition from soft to hard physics: 
at low energies the soft contribution dominates, at high energies the hard and semi-hard 
ones, at intermediate energies (that is where experiments are performed presently) all 
contributions are important. 

Let us consider nucleus-nucleus [AB) scattering. In the Glauber-Gribov approach 
||^, the nucleus-nucleus scattering amplitude is defined by the sum of contributions 
of diagrams, corresponding to multiple elementary scattering processes between parton 
constituents of projectile and target nucleons. These elementary scatterings are exactly 
discussed above, namely the sum of soft, semi-hard, and hard contributions: T2^2 = 
Tsoit + Tscmi + Thard- A corrcspondiug relation holds for the inelastic amplitude T2^x- We 
use the above definition of a cut elementary diagram, which is graphically represented 
by a vertical dashed line, whereas the elastic amplitude is represented by an unbroken 
line: 



2^2. 



This is very handy for treating the nuclear scattering model. We define the model via 
the elastic scattering amplitude Tab~^ab which is assumed to consist of purely parallel 
elementary interactions between partonic constituents, as discussed above. The ampli- 
tude is therefore a sum of terms as the one shown in fig. |5[ One has to be careful about 
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Figure 5: Elastic nucleus-nucleus scattering amplitude, being composed of purely parallel elementary Inter- 
actions between partonic constituents of the nucleons (blobs). 



energy conservation: all the partonic constituents (lines) leaving a nucleon (blob) have 
to share the momentum of the nucleon. So in the explicit formula one has an integra- 
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tion over momentum fractions of the partons, taking care of momentum conservation. 
Having defined elastic scattering, inelastic scattering and particle production is practi- 
cally given, if one employs a quantum mechanically self-consistent picture. Let us now 
consider inelastic scattering: one has of course the same parallel structure, just some 
of the elementary interactions may be inelastic, some elastic. The inelastic amplitude 
being a sum over many terms - Tab^x = J2i "^ab^x ~ ^° squared and summed 
over final states in order to get the inelastic cross section, which provides interference 

terms J2x('^AB^x)i'rAB-*xT ' which can be conveniently expressed in terms of the cut 
and uncut elementary diagrams, as shown in fig. |^. So we are doing nothing more 
than following basic rules of quantum mechanics. Of course a diagram with 3 inelas- 
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Figure 6: Typical interference term contributing to the squared inelastic amplitude. 

tic elementary interactions does not interfere with the one with 300, because the final 
states are different. So it makes sense to define classes K of interference terms (cut 
diagrams) contributing to the same final state, as all diagrams with a certain number of 
inelastic interactions and with fixed momentum fractions of the corresponding partonic 
constituents. One then sums over all terms within each class K, and obtains for the 
inelastic cross section 

where we use the symbolic notation d^b ~ J (Pbo J (P^bA p{bA) J cP^bs p{bB) which means 
integration over impact parameter bo and in addition averaging over nuclear coordinates 
for projectile and target. The variable K is characterized by AB numbers nik represent- 
ing the number of cut elementary diagrams for each possible pair of nucleons and all 
the momentum fractions a;+ and of all these elementary interactions (so if is a partly 
discrete and partly continuous variable, and J2 meant to represent J2 /)• This is the 
really new and very important feature of our approach: we keep explicitly the depen- 
dence on the longitudinal momenta, assuring energy conservation at any level of our 
calculation. 

The calculation of fl actually very difficult and technical, but it can be done and we 
refer the interested reader to the literature ||^] . 

The quantity (if) can now be interpreted as the probability to produce a configu- 
ration K at given s and b. So we have a solid basis for applying Monte Carlo techniques: 
one generates configurations K according to the probability distribution fl and one may 
then calculate mean values of observables by averaging Monte Carlo samples. The prob- 
lem is that fl represents a very high dimensional probability distribution, and it is not 
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obvious how to deal with it. We decided to develop powerful Markov chain techniques 
1^] in order to avoid to introduce additional approximations. 



3 Markov Chain Techniques 

In order to generate K according to the given distribution ^{K), defined earlier, we 
construct a Markov chain 

such that the final configurations X^*"""^) are distributed according to the probability 
distribution fi {K) , if possible for a t^ax not too large! To obtain a new configuration 
= L from a given configuration = K. We use Metropolis' Ansatz for the tran- 
sition probability L) as a product of a proposition matrix w{K, L) and an acceptance 
matrix u{K,L). The detailed balance condition - which assures the convergence of the 
chain - is automatically fulfilled if u{K, L) is defined as 

We are free to choose w{K, L), but of course, for practical reasons, we want to minimize 
the autocorrelation time, which requires a careful definition of w. An efficient procedure 
requires u{K,L) to be not too small (to avoid too many rejections), so an ideal choice 
would be w {K, L) = Q {L). This is of course not possible, but we choose w{K, L) to be a 
"reasonable" approximation to Q,{L) UK and L are reasonably close, otherwise w should 
be zero. So we define 

^ ' [ otherwise 

where d{K, L) is an integer quantity representing a distance between two configurations 
(the maximum number of elementary interactions being different!^), and where Oo has 
a simple structure, just being a product of terms po representing one single elementary 
interaction. So one proposes only new configurations being "close" to the old ones. The 
above definition of w{K, L) may be realized by the following algorithm: 



choose randomly a particular elementary interaction (say the /i*'' interaction of the 
fc**^ nucleon-nucleon pair) 



• propose a new configuration L, which is obtained from the old one K by removing 
the interaction of the fc*^ nucleon-nucleon pair, and replacing this by a new one 
according to the distribution po • 

This proposal is the accepted with a probability u{K, L). One should note that proposing 
a configuration according to some "approximation" fig of il is fully compensated by the 
acceptance procedure, so it is an exact numerical solution of the problem, whatever be 
the precise definition of fio . 

This procedure works extremely well. We performed many test for situations where 
conventional techniques work as well, and we find excellent agreement by using 66.7 fcmax 
iterations, where fcmax is an upper limit estimate of the number of nucleon-nucleon 
interactions. The Markov-chain method is perfect for our purposes, because we have 
fast convergence due to the fact that f2o is not too different from Q., on the other hand 
one cannot use just Oq to obtain an approximate solution, because here we introduce 
an substantial error, which reaches for example already on the nucleon-nucleon level 
about 100 %. 
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4 Summary 



We provide a new formulation of the multiple scattering mechanism In nucleus-nucleus 
scattering, where the basic guideline is theoretical consistency. We avoid in this way 
many of the problems encountered in present day models. We also introduce the neces- 
sary numerical techniques to apply the formalism In order to perform practical calcula- 
tions. 

This work has been funded in part by the 1N2P3/CNRS (PICS 580) and the Russian 
Foundation of Basic Researches (RFBR-98-02-22024}. 
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